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We investigate means to describe the non-local properties 
of quantum systems and to test if two quantum systems are 
locally equivalent. For this we consider quantum systems that 
consist of several subsystems, especially multiple qubits, i.e., 
systems consisting of subsystems of dimension two. We com- 
pute invariant polynomials, i.e., polynomial functions of the 
entries of the density operator which are invariant under local 
unitary operations. 

As an example, we consider a system of two qubits. We 
compute the Molien series for the corresponding representa- 
tion which gives information about the number of linearly in- 
dependent invariants. Furthermore, we present a set of poly- 
nomials which generate all invariants (at least) up to degree 
23. Finally, the use of invariants to check whether two density 
operators are locally equivalent is demonstrated. 



I. INTRODUCTION 

Non-locality is one of the astonishing phenomena in 
quantum mechanics. Well-known examples are EPR 
pairs 1^ and the GHZ state ^ . States of quantum codes 
contradict local realism, too ||^. One common feature 
of these states is that the non-local properties do not 
change under local transformations, i.e., unitary oper- 
ations acting independently on each of the subsystems. 
Thus, any function invariant under local unitary trans- 
formations (LUT) can be used to describe these non-local 
properties p^jl^ . Here we study polynomials that are 
invariant under local unitary transformations. Among 
these, there are e. g. the coefficients of the characteristic 
polynomial of a density operator or of the reduced den- 
sity operators. The paper extends the work of Rains p^ ] 
in which he showed how, in principle, all polynomial in- 
variants can be computed. We present further reductions 
of complexity which make the computation of polynomial 
invariants feasible at least for small systems. 

The paper is organized as follows: In Section || we con- 
sider the linear action and the action by conjugation of 
matrix groups on polynomials. Then we establish a con- 
nection between invariant polynomials and the algebra 
of matrices commuting with all elements of the group. 
A physical interpretation of the invariant polynomials 
is given by relating them to some observables. Mainly 



classical results for these algebras are recalled in See- 
In Section IV a method to construct a vector 
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space basis of these algebras is presented for the case of 
two-dimensional subsystems. Furthermore, we present 
results that imply a further reduction of the complexity 
to compute all invariants. The special situation of pure 
states and subspaces is considered in Section In Sec- 
tion VI we compute the Molien series and a set of invari- 
ants for a two qubit system. We conclude in Section VII 
with examples for the application of these invariants. 



II. POLYNOMIAL INVARIANTS 

A. Operation on polynomials 

The group GL{n,F) of invertible nxn matrices over 
the field F operates on the polynomials p{xi, . . . , Xn) G 
F[xi, . . . , Xn] in the following manner: 



p^{xi, ...,x„) :=p(a;i, . . . ,a;„) 

where (xi, . . . ,5„) := {xi, . . . ,a:„)9 



(1) 



i.e., each variable is substituted by the linear combina- 
tion obtained by multiplying the vector of variables by 
the group element g £ GL(n,F). 

On nxn matrices the group GL{n, F) acts by conjuga- 
tion. Hence polynomials f{pij) = /(pii, ■ • ■ , Pnn) in the 
entries pij oi a n x n density operator p, are then acted 
upon by conjugation, i.e.. 



/^(Py ) = /(Py ) where p := p9 = 



P-9- 



(2) 



Given a subgroup G < GL(n, F), we are interested in 
polynomials that are fixed by all elements of G under 
the action defined by either (|l|) or (||). These invariant 
polynomials (just called invariants) form an algebra over 
the field F since any linear combination and any product 
of invariants is invariant under the action of the group, 
too. It is sufficient to study homogeneous polynomials 
as each homogeneous polynomial of degree k remains ho- 
mogeneous of the same degree under the operation of G 
and every polynomial can be decomposed additively into 
its homogeneous components. 

For the class of so-called reductive groups (e. g. finite 
groups, unitary groups) the invariant ring is finitely gen- 
erated (cf. 0), i. e., every invariant can be expressed in 
terms of some algebra generators. These so-called fun- 
damental invariants can be chosen to be homogeneous 
polynomials of small degree. Under this assertion the 
task is to find a system of fundamental invariants such 
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that any other invariant can be expressed as a polyno- 
mial of these. In what follows we focus on this task for 
invariants under the action of tensor products of unitary 
groups on density operators by conjugation given by (H). 



is very difficult to perform. In Section III we will 



present a method to calculate the matrices F directly 
without integration. 



B. Invariant matrices 



Physical interpretation of the polynomial 
invariants 



Instead of studying the invariant polynomials directly, 
we use the relation between homogeneous polynomials in 
the entries of a density operator p and constant matrices. 

Lemma 1. For every homogeneous polynomial / of 
degree k in the entries of the density operator p there 
exists a matrix F such that 



nP^3) = fF{p^,) = tr (F • p®'^) 



(3) 



Proof. This follows directly from the fact that the ma- 
trix p*^*^ contains all monomials in the variables pij of 
degree k. 

□ 

Next, we characterize invariant polynomials in terms 
of the corresponding matrices. 

Theorem 2. A homogeneous polynomial / of degree 
k in the entries of the density operator p is invariant 
under the operation of a compact group G < GL(n,F) 
by conjugation if and only if f — fp for a matrix F that 
is invariant under conjugation by {g~^)'^'' (equivalently, 
iff the matrix F commutes with 5®'^) for all g € G. 

Proof. Conjugation of p by g corresponds to conjuga- 
tion of F by {g~^)^'' as shown by the following calcula- 
tion: 



If {{g ^ ■ p-g)ij) 



tr(F.(g-i-p.5)®'=) 
tr(F-(g-i)®'=-p®'=-5®'=; 
trig®'' - F ■{g-^)®'' -p®^) 

^ V- ' 

F 

fpiPij)- 



If F commutes with g®'', then the equality of F and F 
implies fp {{g~^pg)ij) = fF{pij)- 

If on the other hand / is invariant under the oper- 
ation of G, for any matrix F with f = fp we have 
/f {{g~^pg)%j) = fFiPij)- Since the group G is com- 
pact, we can average over the group (cf. [^o)) and obtain 
the matrix 



F= I (g®'' ■ F ■ [g-T") dpGig). 
'geG 



(4) 



By construction, F is invariant under conjugation by 
{g"'^)®'' and furthermore / = f-p. 

□ 

Using this theorem and Lemma 1, we are in principle 
able to compute invariants of the group G starting from 
any matrix F and computing a matrix F that commutes 
with g®'' for all g ^ G. But in practice, the integration 



Although we do not have full insight in the physical 
interpretation of the polynomial invariants yet, we will 
relate them to some observables. 

Recall from Eq. (^) that all polynomial invariants of 
degree k can be written as f{pij)F = tr (F • p*^*^). From 
F we construct two Hermitian operators 

Ml := F + F'' and M2 -.^ i F ~ i F'' . 

Both Ml and M2 commute with g®'' for all g = Ui (gi 
... Un e U{n)®'^ since F (and thus F'^) commutes 
with g®'^ (cf. Theorem 2). Hence the (real) mean values 

(Afi) tr {Ah ■ P®'=) and (Ah) tr {Ah ■ p®"") 

are also invariant under local unitary transformation. In 
principle they can be obtained by joint measurements of 
k copies of the quantum system with density operator p. 



III. INVARIANT ALGEBRAS 

In order to compute all homogeneous invariants of de- 
gree k it is sufficient to know the algebra of matrices 
that commute with g®'' for all g E G. Such algebras 
have been studied e. g. by Braucr for many classes of 
groups (cf. 1^). For the unitary group and tensor prod- 
ucts of unitary groups, we have the following theorems 
and corollaries: 

Theorem 3. (Brauer) The matrix algebra An.k of ma- 
trices that commute with any matrix J7®'^ for U G U{n) 
is generated by the representation Tn,k'. Sk GL{n'',C) 
of the symmetric group Sk that operates on the tensor 
product space (C")®*^ = Vi ® . . . (X) Vfe by permuting the 
k spaces Vi of dimension n. 

This result extends to tensor products of unitary 
groups: 

Corollary 4- The algebra of matrices that commute 
with any matrix i/f" ... U'^'' for Ui eU{ni) is given 
by the tensor product of the algebras Am^k- 

To obtain the "natural" ordering of the tensor factors, 
we have to conjugate the matrices by a permutation ma- 
trix: 



of matrices that com- 



Corollary 5. The algebra A^^ 

mute with any matrix {Ui ® . . . Un)^'' for Ui E U{n) 
is conjugated to the iV-fold tensor product of the algebra 
An J /c , i.e.. 



A. 



(N) 



a {A„ 



kN 



(r). (5) 
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Here r is the permutation that exchanges the macro- and 
micro-coordinates according to the isomorphism between 
the tensor product spaces 

(V^^f and {V'^^f (where y = C"). 

As a permutation on {1, . . . ,k ■ N}^ r maps ak + b + I to 
67V a + 1 (for a = 0, . . . , iV - 1, 6 = 0, . . . , fc - 1). 

(The reader famihar with the theory of fast Fourier 
transformations will recognize the similarity to the "bit 
reversal permutation" ||].) 

For the special situation of qubit systems, i.e., the 
group U{2), the dimension of the algebra is given by the 
following theorem: 

Theorem 6. The vector space dimension of the algebra 
•A.2,k is given by the Catalan number |^ 

Proof. This result is derived in jl^ from a theorem of 
Weyl 0. □ 

Note that the algebra ^2,fc was defined by the k\ matri- 
ces r2,fc(7r) for TT e Sk- As an algebra, A2,k is generated 
by the image of the generators of Sk, i.e., by only two 
matrices. Theorem 6 states that, even as a vector space, 
less than k\ matrices are sufficient. 



IV. BINARY TREES, PERMUTATIONS, AND 
ALGEBRAS 

A. One qubit 

The mapping (^ from invariant matrices to invariant 
homogeneous polynomials is a vector space homomor- 
phism. Thus, in order to compute all linearly indepen- 
dent homogeneous invariants of degree k for tensor prod- 
ucts of the group J7(2), it is sufficient to consider a vec- 
tor space basis of the algebra A2,k- Such a basis can 
be constructed starting from binary trees with k nodes, 
mapping them to permutations of k letters, and finally 
obtaining matrices via the representation T2^k- The con- 
struction resembles some of the many beautiful combina- 
torial properties of Catalan numbers (cf . ) . 

Let Bk denote a labeled ordered binary tree with k 
nodes, i. e., each node in the tree but the root has a 
father, and each node in the tree has at most one left 
and at most one right son. The labeling of the k nodes 
of the tree with the numbers {f , . . . , A:} is obtained by 
traversing the nodes in the order root, left sub-tree, right 
sub-tree. FIG. |l| shows all distinct binary trees with three 
nodes labeled in that manner. 

A maximal right path in the binary tree Bk is a se- 
quence of nodes (rg , ri , . . . , ) such that each of the 
nodes r.i+i is the right son of the node r^, rg is not the 
right son of any node, and rj has no right son. 




FIG. 1. The five distinct ordered binary trees with three 
nodes. 

Given the set Tl{Bk) of all maximal right paths of a 
binary tree Bk we define a permutation TT{Bk) G Sk by 
the product of cycles 

^{Bk)^ n (rori...r,J. (6) 

{ro,ru...,rj.)eniBk) 

E. g., for the trees of FIG. |^ we get the five permutations 
(1)(2)(3), (1)(2 3), (13)(2), (1 2)(3),and (12 3). 

Let Bk denote the set of all distinct ordered binary 
trees with k nodes labeled in the manner described be- 
fore, and be Vk = i^{Bk) the set of permutations obtained 
by the mapping tt. Using this notation, we can formulate 
the following theorem: 

Theorem 7. The set of matrices Mk '■— T2.k{Vk) = 
{T2.k{T^{B)) : B G Bk} forms a vector space basis of the 
algebra A2,k- 

Proof. For k = and k — 1 the statement is obviously 
true. 

For A: > 1 , we partition the set Bk of binary trees with 
k nodes into k classes Bk.j {j — 0, . . . , k — 1). The class 
Bk,j consists of all trees with j nodes in the left sub-tree 
of the root. The general form of a tree in the class Bk,j 
is shown in FIG. B. 




FIG. 2. General form of a tree in the class Bk,j- The 
sub-trees Bn, Bir, Bri, and Brr might be empty. Further- 
more, for j — (or j — k ~ 1) the left sub-tree Bi (or the 
right sub-tree B,-, resp.) is empty. 

For j < fc ~ 1 in each of the trees B e Bkj, j + 2 
is a right son of 1, and thus tt{B) maps 1 to j -I- 2k. 
For j — k — 1, the root 1 has no right son and thus 
tt{B) fixes 1. To combine these two cases, we identify 
k + 1 and 1. Hence T2^k{T^{B)) maps |ei) to |ej+2) where 
\ei) = \0f'-^ |1) 10)®''"'. This shows that for B G Bkj 

tr{\e,){e,,+2\T2MAB))) =Sj,j,. 

Therefore, for j ^ j' the matrices in the sets T'2.fc(7r(Sfej )) 
and T2.k{'^{Bk,j')) are mutually linearly independent. 

For fixed j, each permutation ■k{B) for a tree B G Bk,j 
with left and right sub-trees Bi and B^ (see FIG. ||) can 
be written in the form 

tt{B) = (1 j + 2) • Tl{Bi) ■ n{Br) 
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where the permutations tt; = 7r(_B/) and tt,. = n{Br) op- 
erate on the sets {2, . . . , j + 1} and -f 2, . . . , fc}, resp. 
The corresponding representations are "shifted" tensor 
products, i. e., 

T2m{AB)) = T2,fe((l J + 2)) • (I2.+1 ® T2^fc_j_iK)) 

■(I2I ®T2j(7r;') ® l2fc-.-l) 

= T2,fc((lj + 2)) 

•(I21 ®T2,,(7r;')®T2,fc_,_i(7r;)) 

where tt; G 5^ and tt^ G Sk-j-i are obtained by rela- 
behng. By the induction hypothesis, the matrices cor- 
responding to the sub-trees Bi and are hnearly in- 
dependent for different trees. Thus we have shown that 
the matrices in Aik are Hnearly independent. It remains 
to show that they form a basis of A2,k- But this follows 
from Theorem 6 together with the fact that there are 
exactly C{k) different ordered binary trees with k nodes 
(cf. p. 389]). 

□ 



B. Multiple qubits 

In order to compute a basis of the algebra ^4^^^^ for an 
N particle system, we define the following representation 
Tf!k- {Sk)^ GL{n^^, C) of the iV-fold direct product 
of the symmetric group Sk'- 

(tTi, . . . , TTat) 1-^ cr • (T„,fc(7ri) ® Tn,k{l^N)) ■ o-"^ 

(where the matrix a is given by (|^)). 

In the special case of N qubits (i. e. n = 2) combining 
Theorem 7 and Corollary 5 we obtain: 

Corollary 8. The set of matrices T2^\'P^) is a vector 

space basis of the algebra -42^,'' . 

So far, we are able to compute a vector space basis for 
A2^i} as follows: 

• Generate the set Bk of all different binary trees. 

• Generate the set of permutations Vk obtained by 
construction (^. 

• For each iV-tuple of permutations (tti, . . . , ttjv) ap- 
ply the representation Tj^^ i.e., compute the ten- 
sor product of the representations T2^ki'^v)- 

Instead of computing a matrix for each of the (fc!)^ tu- 
ples of permutations in (Sk)^ , it is sufficient to consider 
only the C{k)^ permutations, implying a complexity re- 
duction from 0{k''^) to 0(4*=^). 

Using equation (||), we get a set of polynomials in- 
variant under local transformations spanning the vector 
space of homogeneous polynomial invariants of degree k. 

For any A^-tuplc vr = (tti, . . . , ttjv) G (Sk)^ of permu- 
tations we obtain a homogeneous invariant of degree k 
given by 



Clearly, there exist relations between the invariant 
polynomials obtained from the tuple of permutations 
(tti, . . . , ttat) e (Sk)^ (also with varying k). Some of 
these relations can be expressed in terms of the permu- 
tations. This allows a further reduction of the number 
of tuples of permutations that have to be considered to 
compute all invariants. 

Theorem 9. If (tti, . . . , npf) and {Tr[, . . . , tt^) are "si- 
multaneously" conjugated, i.e., there exists r G Sk 
such that ttI = t^-^tTi,t for all ly G {!,..., A^}, then 

/7ri,...,7rjv (Pij) — fir[,...,Tr'^iPij)- 

Proof. Simultaneous conjugation of the permutations 
t:^ by T corresponds to permuting the tensor factors 
in p®'^ by r, keeping p^'' fixed as a whole. Thus 
/7ri,...,7rjv(Pii) does not change. 

□ 

Next we give a condition on the permutations when an 
invariant can be written as a product of invariants. 

Theorem 10. If the subgroup H < Sk generated by 
TTi, . . . , TT^r is not trausitlve, then the homogeneous invari- 
ant /7ri,...,7rjv (Pij) of degree fc is a product of invariants of 
smaller degree. 

Proof. If H is not transitive, it defines a non-trivial 
partition of the set {1, . . . , A:} into orbits. By simultane- 
ous conjugation using Theorem 9, we can assume that the 
partition {1, . . . , fci} and {fci -|- 1, . . . , /c} respects the or- 
bits. Thus, for V — 1^ . . . , N each permutation tt^ can be 
written as a product -k^, — tt^ - tt" with 7r[, G Ski operating 
on {1, . . . , fci} and tt" G operating on {fci + 1, . . . , fc} 
(fci + ^2 — k). Furthermore, 

/.,,...,.„(p.,)=tr(riy(7ri,...,7r^).p®'=) 

= tr((TW(^'i, ...M® . . .yk)) ■ P^') 

= tr(ri2(7ri,...,.5v)-p^'=0 

.tr(rW(7ri',...,7r^).p«'=^) 

= fi:[,...,TT'^{Pii) • fi:'(,...y^{Pij)- 

□ 

For the case of two qubits. Table | shows the reduction 
of the number of pairs of permutations to be considered 
using the construction of permutations from binary trees. 
Theorem 9, Theorem 10, both Theorem 9 and 10, and 
finally the combination of Theorems 9, 10, and 11. 

V. PURE STATES AND SUBSPACES 

The technique to compute polynomials invariant under 
the action of tensor products of unitary groups does not 
only apply to density operators of mixed states, but also 
to subspaces and pure states. To study non-local prop- 
erties of subspaces with basis |V'i) (e.g. quantum error- 
correcting codes), one can use the invariants of the cor- 
responding projection operator P = Pure 
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Th. 9 


Th. 10 


Th. 9 
Th. 10 


Th. 9 
Th. 10 
Th. 11 


1 


1 


1 


1 


1 


1 


1 
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36 


25 


10 


15 


6 


3 


4 


(J 1 u 


1 Qfi 


ou 


Q7 




1 n 


5 


14 400 


1764 


114 


733 


60 


22 


6 


518 400 


17424 


496 


6147 


291 


100 


7 


25 401 600 


184 041 


2142 


55 541 


1310 


361 


8 


1 625 702 400 


2 044 900 


10 758 


530 773 


6 975 


1717 



■ j(2) ■ . . . ■ Pj(7r(k)) J(k) . (8) 



TABLE I. Number of pairs (7ri,7r2) G {Sk)^ to be con- 
sidered for the construction of invariants using the different 
theorems. 



states \(f)) can be considered as one-dimensional subspace 
with projection operator P = |0) or equivalently as 
a mixed state with density operator p = 

In that situation, we have the additional relation — 
P (resp. = p) which can be used for a further reduc- 
tion of the number of permutations to be considered. The 
following theorem is quoted from [p^ , adding an explicit 
proof. 

Theorem 11. Let P be a projection operator. If for 
(tti, . . . , ttat) G [Sk)^ there exist different numbers I and 
TO such that for each permutation tTj, we have t:^{1) = to, 
then the invariant /7ri,...,7r„(fy) = f .y^iPi]) where 
the permutations tt^ 6 Sk-i are obtained by identifying 
the points I and to followed by a relabeling. 

Proof. By Theorem 9, we can assume without loss of 
generality that / = 1 and to = 2, i.e., 7r^(l) = 2 for all 
permutations 7r,y. We will show that in the summation 
there are entries of P^ which can be replaced by those 
of P. The entries of P®*^ are of the form 

Here the indices i^^^ are A^-tuples ■ ■ ■ ,ijv^) (^'^'' ^ 
system with N particles). The subscript of li'"^ corre- 
sponds to the i/-th particle, whereas the superscript cor- 
responds to the p-th copy of the whole system. Left- 
multiplication by P'^^''(7ri, ... ,772) permutes the rows of 
p®fc yigidijig t^jie matrix M with entries 



Mai 



(i(i),...,i('=)),(j<^',---,i<'=') 



Pi(7r(l)) ■ . . . ■ P^(^(k)) J(k) . 



The operation of tt = (tti, . . . , ttjv) on the indices is given 



by ^('^(''^ 
have 



= Pj(2) j(l) • Pj(x(2)) .,-(2) • 

Now taking the trace results in 

,■(1) ,-(fc) 



Since 7rjy(l) = 2 we 



• Pi(^{k))j(k) . 



Considering the summation over j^^^ separately and us- 
ing P^ = P, we get 



j(,r(2))^j-(2)Pj(2) — Pj(^(2}) . 



(9) 



Combining equations (||) and (||) yields 

E ^(jd),. ..,»('')), (id), ...jC')) 
,...,i('=) 

= ^ P,(,r(2)) .,-(l) • Pj(^(3)) ,j(3) • . . . • Pj[^(k)) j{k) . 

j"(l),j(3),...,i('=) 

Now the result follows immediately if we identify the 
points 1 and 2. 

□ 



VI. THE INVARIANT RING OF A TWO QUBIT 
SYSTEM 

To illustrate the results, we consider the smallest non- 
trivial example, a system of two qubits. Using our algo- 
rithm, we are able to compute homogeneous invariants 
for each degree. As stated before, the homogeneous in- 
variants of fixed degree form a vector space. Therefore it 
is sufficient to compute a basis for that vector space, e. g., 
a maximal linearly independent set of homogeneous in- 
variants. In order to know how many invariants we need, 
we address the problem of determining the dimension dk 
of the vector space of invariants of degree k. Information 
about these dimensions is encoded in a formal power se- 
ries, the Molien series (cf. p^) 

P(z) :=^d,z'=eZ[[z]]. 

fc>0 

In case of a finitely generated algebra the Molien series 
turns out to be a rational function (see e. g. Thus it 
can be expressed in a closed form with a finite number 
of terms. In principle, the Molien series for the linear 
action (|l|) of a compact group G on polynomials can be 
computed by means of the following averaging formula: 



P(z) 



geG det(id 



(10) 



where pc is the normalized Haar measure of G. 

In this paper, however, we are concerned with the ac- 
tion of G = SU{2) X SU{2) on matrices by conjugation 
(^, for which the formula ( [To| ) does not apply directly. 
Since the operation on p given by 

p ^ (C/i (g) C/2) • p • [Ui ® ?72)t for C/i, C/2 e SU{2) 
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is linear, we write p as a vector p of length r? and obtain 
the representation 

where Ui denotes the complex conjugate of the matrix 

The integral (|l^) is simplified by means of the integral 
formula of Weyl (cf. §26.2]) which allows to perform 
an integration over the whole group in two steps. The 
first step involves an integration over a maximal torus 
T of the group G and the second an integration on the 
residue classes G/T. Next, the integral is transformed 
into a complex path integral which can be solved by the 
theorem of residues. Finally, we end up with the Molien 
series 

p z^O - 2* - z'' + + 2z5 + 2z-4 - z3 - z^ + 1 

^^-^ ~ (z - l)10(z + l)6(z2 + l)2(z2 + Z + 1)3 

= 1 + z + 4z2 + 6z3 + 16z'' + 23z5 + 52z'' + llz' + 150z** 
+ 224z^ + 396z^° + 583z" + 964z^2 _,_ 13952I3 
+ 2\mz^^ + SlOOz^^ + 4639zi'' + 6466zi'' + 9344zi* 
+ 12785z^^ + 17936z2f' + 24121z2i + 33008z22 
+ 43674z23 +0{z24). 

The information about the dimensions of the vector 
spaces can now be used to compute all invariants degree 
by degree. Having computed an algebra basis for all in- 
variants of degree less than fc, homogeneous invariants of 
degree fc are obtained by multiplying invariants of smaller 
degrees that sum up to fc. By computing the vector space 
dimension of these invariants and comparing it to the di- 
mension given by the Molien series, we know how many 
linearly independent invariants are missing. Next, these 
missing invariants are constructed from pairs of permu- 
tations "Ku G S'fc. It is sufficient to draw randomly from 
the pairs of permutation remaining after Theorem 9 and 
Theorem 10 until the vector space dimension given by 
the Molien series. 

Using the computer algebra system Magma [|j we 
found 21 invariants corresponding to the permutations 
shown in Table ||. Furthermore, we were able to show 
that these invariants generate all invariants up to degree 
23. We conjecture that any polynomial invariant of a 
two-qubit system can be expressed in terms of these 21 
invariants. It should be noted that there are 43 674 lin- 
early independent invariants of degree 23 and that the 
invariants of degree 9 have more than 24 000 terms. 

The Molien series provides also information about the 
maximal number of algebraically independent invariants. 
This number is given by the order of the pole z = 1 of the 
Molien series [|l^. For the two-qubit system, there are 
10 algebraically independent invariants. Thus any 11 in- 
variants fulfill a polynomial equation, i. e., given numeri- 
cal values for 10 algebraically independent invariants, the 
values of the remaining invariants are some roots of poly- 
nomials. But these values are not unique since none of 
the 21 invariants is a polynomial function of the others. 





degree 




7r2 




^ terms 


1 


id 


id 


4 


2 


id 


(12) 


10 


2 


(12) 


id 


10 


2 


(12) 


(12) 


10 


3 


(12 3) 


(12) 


52 


3 


(12 3) 


(12 3) 


24 


4 


(12 3 4) 


(13) 


110 


4 


(12 3 4) 


(12 3) 


144 


4 


(12 3 4) 


(1234) 


70 


4 


(12 3 4) 


(12)(34) 


98 


5 


(12 3) (4 5) 


(12345) 


456 


6 


(123456) 


(1235) 


1334 


6 


(12 3) (4 5) 


(12456) 


1586 


6 


(12 3) (4 5) 


(123456) 


1542 


6 


(12 3) (4 5) 


(1234)(56) 


1464 


7 


(1234)(567) 


(124567) 


4156 


7 


(12 3 4) (5 6 7) 


(12 6 7)(3 5) 


4 576 


8 


(12345)(678) 


(1235678) 


10 414 


8 


(12345)(678) 


(12378)(46) 


11340 


9 


(12345)(678) 


(1236789)(45) 


24 780 


9 


(12345)(6789) 


(1235678) 


24168 



TABLE II. The invariants .f-ni,TT2{pij) corresponding to 
these permutations (listed together with the degree and the 
number of terms of the invariants) generate the polynomial 
invariants of a two qubit system (at least) up to degree 23. 



VII. EXAMPLES 

A. Characteristic polynomials 

As stated in the introduction, the coefficients of the 
characteristic polynomial Xp(X) of a density operator p 
and of the reduced density operators are invariant under 
local unitary transformation. They can be expressed in 
terms of the invariants /Tri,7r2 presented in Table |l| as 
follows: 

Xp(X) = - fidAdX^ + (^-fid,id - 2),(l 2)^ 

+ i^^-^ffdAd + -fid,idf(12),(l 2) - -/(I 2 3), (1 2 3)^ -'^^ 

+ '^fid,id - jfidAdfil 2), (12) + -/id,w/(1 23),(1 23) 

+ ^f2 1. 

+ g/(12),(12) - J/(1234),(1234)- 

Here, the coefficient of in Xp is a linear invariant 
polynomial that equals the negative trace of p, and the 
constant coefficient of xp is an invariant of degree four 
that equals the determinant of p. 

For the characteristic polynomials of the reduced den- 
sity operators p^^^ := tr2(p) and p^^^ := tri(p) we obtain 

Xp(i) i^) = X'^ ~ fid.idX + 2 ifid,id ^ f(12),id) ^'^'^ 

Xp(2) {X) = - fid,idX + - {ffd,id - /id,(i2)) ■ 
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B. Test for local equivalence 



ACKNOWLEDGMENTS 



Consider the following density operators: 



Pi 



and 



P2 




-16-150i -10 + 8j 

226 13 

13 234 

10 + 8i 16- 150 i 




1 
800 



f 214 

-16 - 150 i 
10 - 8i 
V -3 



-16 + 150i 
226 
-13 

-10 + 8i 



10 + 8i -3 

-13 -10 -8i 

234 16 - 150 i 

16 + 150 i 126 



These density operators are (globaUy) conjugated to each 
other, i. e., p2 = UpiW where U S C/(4). Thus they have 
the same eigenvalues. Furthermore, the reduced density 
operators are the same. Partial transposition of pi and 
P2 yields positive operators pj' and p^* (z = 1,2), thus 
pi and p2 are both separable . Again, the eigenvalues 
and the reduced density operators corresponding to pf' 
and P2^ are the same. Also, all polynomial invariants up 
to degree five evaluate to the same number for pi and 
P2. Hence one might expect that pi and p2 are locally 
equivalent. 

But we have the following: The values of the invariants 
of degree six corresponding to tti — (12 3)(4 5), n2 — 
(12456) and tti = (123)(45), tts = (123456) differ 
showing that pi and p2 are not locally equivalent. Note 
that after pre-computation of the invariants this can be 
decided just by evaluating the invariants. 



VIII. CONCLUSION 

We have established a connection between local poly- 
nomial invariants of quantum systems and permutations. 
The local invariants of an N particle system can be com- 
puted directly from A^-tuples of permutations. Theo- 
rems 9 and 10 allow a reduction of the number of permu- 
tations to be considered. This makes the computation of 
invariants feasible. 

In the particular case of pure states and subspaces. 
Theorem 11 allows a further reduction of complexity. All 
these theorems apply for subsystems of any dimension. 

In the special case of qubits systems, we are able to 
construct a vector space basis for the matrices directly in 
terms of permutations derived from binary trees. It has 
to be investigated if there are similar constructions for 
subsystems of dimension n > 2. 

Since our methods are not restricted to density opera- 
tors, they can also be used to study non-local properties 
of unitary transformations. For example, it can be tested 
whether two quantum circuits are equivalent with respect 
to conjugation by single qubit gates. 
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